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Abstract. In this paper we model pedestrian flows evacuating a narrow corridor 

f^ ' through an exit by a one-dimensional hyperbolic conservation law with a non-local 

^Nj , constraint. Existence and stability results for the Cauchy problem with Lipschitz 

» I ' constraint are achieved by a procedure that combines the wave-front tracking al- 

f~^ l' gorithm with the operator splitting method. The Riemann problem with piecewise 

■^C . constant constraint is discussed in details, stressing the possible lack of uniqueness, 

self-similarity and Ljo^-continuity. One explicit example of application is provided. 

m 
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c^ ■ 1. Introduction 

The theory for constrained conservation laws was introduced by Colombo and 
Goatin in Ref. \8\. Their results are of interest in many real-life applications, such 
as vehicular traffic, Refs. [11], [20], pedestrian flows, Refs. [1], [13], telecommuni- 
cations, supply-chains, etc. 
QQ ■ For pedestrians, constraints are usually caused by a direct capacity reduction 

CN , (door or obstacle) and are of fundamental importance in the calculation of evac- 

uation times. The first macroscopic model for pedestrian evacuations able to 



> 

(N 



^ 



O 



^r . reproduce the fall in the efficiency of an exit when a high density of pedestrians 

clogs it, was the CR model proposed in Ref. pl3] and developed in Refs. [Q], [10], 

|16) . |28j : see also Ref. [29j . There, the maximal outflow allowed through the exit 

is assumed to be a piecewise constant function of the density at the exit, and takes 

_ two distinct values, one related to the "standard" case, when the density is less 

r> ■ than an assigned threshold, and one related to the case with "panic", when the 

C^ , density is greater than the threshold. As a result, the fall in the efficiency of the 

exit has a non-realistic behavior since it is instantaneous when the panic reaches 
the exit. A more realistic model should reproduce a more gradual decay in the 
efficiency of the exit as the pedestrians accumulate close to it, see Refs. |27j, [31j . 
Moreover, according to the CR model, once the efficiency of the exit falls down, it 
remains constant until the very last pedestrian is evacuated. On the contrary, in 
real life the efficiency of the exit gradually increases as the number of the remaining 
pedestrians to be evacuated becomes smaller and smaller. 

To avoid these drawbacks of the CR model, in this paper we generalize the 
results proved in Ref. [8] and study the Cauchy problem for a one-dimensional 
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hyperbolic conservation law with non-local constraint of the form 

dtp + d^f{p) = (t,x)eM+xM (la) 

f{p{t,0±))<p(l w{x) p{t,x) dxj teR+ (lb) 

p{0,x)=po{x) xGM. (Ic) 

Above, p = p{t, x) € [0, R] is the (mean) density at time t G M+ of pedestrians 
moving along the corridor parameterized by the coordinate x G M_. Then, R € M+ 
is the maximal density, / : [0,R] -^ M is the pedestrian flow with pedestrians 
moving in the direction of increasing x, p : M+ -^ M+ prescribes the maximal flow 
allowed through an exit placed in x = as a function of the weighted average 
density of pedestrians in a left neighborhood of the exit, w : M_ -^ M+ is the 
weight function used for the average density and po ^ I^ ^^ [0, -R] is the initial 
(mean) density. Finally, p{t,0—) denotes the left measure theoretic trace along 
the constraint implicitly defined by 
-^ r+oo rO 



lim- / / \p{t,x) - p{t,0-)\ (l){t,x) dx dt = 

eJ-O e Jo J-e 

for all (p € C^(M^;M). The right measure theoretic trace, p{t,0+), is defined 
analogously. 

Observe that if w is regular enough and linix'^-oo u!{x) = 0, then the quantity 



e(t)= / wix) p{t,x)dx (2) 

is the solution of the following Cauchy problem for an ordinary differential equation 
at)= [ w{x) [f{p(t,x))-f{p{t,0-))] dx, m= [ w(x)po{x)dx. 



In real life, when a very high density of pedestrians accumulate near the exit, 
the outgoing flow can be very small, but remains strictly positive. For this reason, 
the efficiency of the exit p is assumed to be always strictly positive. We assume 
that the weight w is an increasing function with compact support because the 
efficiency of the exit is more affected by the closest high densities, while it does 
not take into account "far" densities. In summary, we assume that: 
(F) / € Lip{[0,R]; [0,+oo[), /(O) = = f{R) and there exists p G ]0,R[ such 
that f'{p) {p- p) > for a.e. p G [0, R\. 
(W) w G L°°(R_;M_|_) is an increasing map, ||'U^|Ili(r_r ) = 1 ^^^ there exists 

i^o > such that w{x) =0 for any x < — i^. 
(PO) p takes values in ]0, f{p)] and is a non-increasing map. 

Observe that f{p) < f{p) for any p i^ p and ^(t) E [0, i?], see Fig. [TJ In the 
present work we do not take into account the presence of a panic regime since 
in (F) we assume that the fundamental diagram \p i-^ (P)/(/o))] is bell-shaped. 
Indeed, the CR model introduces a flux that results from the juxtaposition of 
two bell-shaped sub-fluxes corresponding to the two regimes quiet-panic and, 
therefore, does not satisfy the condition (F). The latter assumption (PO) is the 
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Figure 1. Examples of functions satisfying conditions (F), (W), 
(PO) and (P2). 



minimal requirement for ([T]) to be meaningful in the sense of distributions, see 
Definition [TJ While existence for the Riemann problem is proved for piecewise 
constant p, see (P2) in Sec. IH we strengthen the assumption on p to a Lips- 
chitz continuity hypothesis when dealing with the Cauchy problem, see (PI) and 
Theorem [1] in Sec[2j 

We give the definition of solution for problem with nonlocal constraint ([TJ by 
extending the definition of entropy weak solution for a constrained Cauchy problem 
of the form 

(3a) 
(3b) 
(3c) 

Definition 1. Assume conditions (F), (W), (PO). A map p G L°°(M+ x 

M; [0, i?])nC^(M+;Lj-^^(R; [0,i?])) is an entropy weak solution to ([TJ if there exists 
q £ L°°(M_|_; [0, /(p)]) such that the following conditions hold: 

ioom2.i 



dtp + dJip) = 


(t,x) € 


fip{t,0±))<q{t) 


t e 


p{0,x) = po{x) 


X € 



(1) For every test function (p € C^(M , M^-) and for every k € [0,R] 



p - k\dt4> + sign(/9 - k) {f{p) - f{k)) 9^0] dx dt 



+ 2 



/(A;)(/<(t,0)dt 



/(P). 
Po{x) — k\ 0(0, x) dx > , 



and 



f{p{t,0±))<q{t) fora.e.te 
(2) In addition q is linked to p by the relation 



Q{t)=p( f 

\Jr. 



w{x) p{t,x) dx ) for a.e. t € 



(4a) 
(4b) 
(4c) 

(4d) 
(5) 



If q is given a priori, then ((3|) is the definition of entropy weak solution to 
problem ([3]). This item is precisely the Definition 2.1 in Ref. [1], which is a minor 
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generalization of the original Definition 3.2 introduced in Ref. [8]. We refer to 
Proposition 2.6 in Ref. [Ij for a series of equivalent formulations of conditions (JH). 

The lines (|4a|) and (|lc|) originate from the classical Kruzkov Definition 1 in 
Ref. [2^, in the case of the Cauchy problem with no constraints. Lines ()4bp 
and (|4dp account for the constraint. Let us stress that both left and right traces 
at X = of an entropy weak solution exist (see, e.g., Theorem 2.2 in Ref. [1\ which 
is a reformulation of the results of Refs. |26] . |33|). 

Our main result is well-posedness for the nonlocal problem ([T]), see Theorem [TJ 
We show that under the Lipschitz continuity assumption on p there exists a semi- 
group (5j)t>o on L°°(R; [0, i?]) such that p{t,-) = St{po) is the unique solution 
of (dl), and depends continuously on t and po with respect to the Lj^j^-distance. 

The uniqueness result for ([T|) is a consequence of a stability estimate for the 
problem with local constraint ([3]) with respect to the Lj-^^-distance, of the rela- 
tion ([5|) and of the Gronwall inequality. 

The existence result for ([1]) is achieved through an operator splitting method, 
Refs. [6], [7], [11], [15], [18], coupled with the wave-front tracking algorithm, 
Refs. [8], [T7]. This procedure is chosen for two reasons. First, wave-front track- 
ing schemes are able to operate also in the case with panic, when nonclassical 
shocks away from the constraint have to be taken into account, see Ref. [3]. Sec- 
ond, the operator splitting procedure allows us to approximate our problem with 
a problem of type ([3]), namely with a "frozen" constraint. This greatly simplify 
our work because it avoids the difficulties coming from the Riemann solver for the 
nonlocally constrained problem. Indeed, we underline the fact that, differently 
from the constrained Cauchy problems studied in Refs. [I], [8], the maximal flow 
at the constraint for ([1]) depends on the solution itself and, in general, it is an un- 
known variable of the problem. As soon as p is discretized, i.e. p is approximated 
by piecewise constant functions, the solution of the corresponding nonlocally con- 
strained Riemann problem may fail to be unique, Ljjjj,-continuous, consistent and 
self-similar, as we will show in Sec. [H Furthermore, using wave-front tracking for 
a nonlocal problem is quite delicate because one cannot merely juxtapose local 
solutions of Riemann problems, see Remark [T] 

The use of wave-front tracking approximation requires the BV functional set- 
ting. As already observed in Ref. [5], the constraint may cause sharp increases 
in the total variation TV(p) of the solution. To overcome this difficulty, as in 
Refs. [5j, [8^, LUj, [32j . we rather estimate the total variation of ^ o p, where 



^ip) = sign(/, - p) ifip) - f{p)) = / fir) dr . (6) 

Jp 

We stress that ^ is one-to-one, but possibly singular at p = p. Indeed, if p is in 
BV, then also ^ o p is in BV, while the reverse implication does not hold true. 
Therefore we introduce the set 

V={peL^{R;[0,R]) : ^(p) gBV(R;M)} . (7) 

In the first step of our construction, we exploit the finite speed of propagation 
property for the conservation law (jlap to define St on the domain D by coupling 
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the operator splitting method with the wave-front tracking algorithm, see Sec [3l 
Then we extend St to L°° by a density argument, see the second part of the proof 
of Theorem [TJ 

The paper is organized as follows. Sec. [2] and [3] are devoted to the constrained 
Cauchy problem. In Seed] we study ([U with a Riemann initial datum. In Sec. Owe 
apply the model ([T]) to describe the evacuation of a corridor through an exit placed 
at X = 0. All the technical proofs are in Sec. O Conclusions and perspectives are 
outlined in Sec. [71 



2. The Cauchy problem with nonlocal constraint 

In this section we consider the Cauchy problem ([1]) under the hypotheses (F), 
(W) and the following assumption on p: 

(PI) p belongs to Lip ([0, R] ; ]0, f{p)]) and it is a non-increasing map. 

Let us start with the basic properties of entropy weak solutions to ([1]). 

Proposition 1. Let [t >—?■ p{t)] be an entropy weak solution of ([TJ in the sense of 
Definition [TJ Then 

(1) It is also a weak solution of ([Taj), ([Tc]l . 

(2) Any discontinuity satisfies the Rankine-Hugoniot jump condition. 

(3) Any discontinuity away from the constraint is classical, i.e. satisfies the Lax 
entropy inequalities. 

(4) Nonclassical discontinuities, see Refs. [21], [29], may occur only at the con- 
straint location x = 0, and in this case the flow at x = is the maximal 
flow allowed by the constraint. Namely, if the solution contains a nonclassical 
discontinuity for all times t a I, I open in M+, then for a.e. t in I 



f{p{t,0-))=f{p{t,0+))=pU w{x) p{t,x)dx\ . (8) 

Proof. By taking k = 0, then k = R, in @, we deduce that any entropy weak 
solution to ([!]) is also a weak solution to p^ . ([Tcj). As a consequence, p satisfies the 
Rankine-Hugoniot jump condition and, in particular, f {p{t,0—)) = f {p{t,0+)). 
By taking in Q a test function with support in M_|- x M_, then in M+ x M_|_, we see 
that p is also a classical Kruzkov solution to (|lap . (fTc]l in M-|- x M± and therefore 
the jumps in p located at x ^ satisfy the Lax entropy inequalities. Finally, 
we prove property ([8]). This property was observed at the level of problem ([3|), 
see in particular the description of the "germ" Qp in Ref. [I], but it was not 
explicitly stated in the works Refs. [1], [8] devoted to problem ([3]). For the sake of 
completeness, we give an explicit proof of property ([8]) for problem Q. Consider 
the test function 

(j){t,x)= I Ss{z)dz i,{t) , 



" r+ 


oo 








' 


/ 
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where ■0 £ C^(M;M+) is such that ip{0) = 0, while 6^ is a smooth approximation 
of the Dirac mass centered at 0+, 6q_^_, namely 

5e G C^(M;M+), e G R+, supp((5,) C [0,e], ||5e|lLi(R;R) = 1, ^e ^ C- (9) 
Observe that as e goes to zero 

(/)(0, x) = ^ , dt(l){t,x) = 



oo 

6eiz) dz 

\x\-e 



1p{t) -^ 



Then, if we take k = p and </> as test function in Q, we obtain as e goes to zero 

[^ {p{t, 0+)) - ^ (p(t, 0-))] ^l^{t) dt + 2 [ [f{p) - p im)] ^{t) dt > , 

JR+ 

where ^ is defined by ([2]). For the arbitrariness of tp, we have for a.e. t > 

vl; (p(t, 0+)) - ^ (/>(t, 0-)) + 2 [/(p) - p im)] > . 

Therefore, if for t G / the solution has a nonclassical discontinuity at the constraint 
location x = 0, then by the assumption (F) and the Rankine-Hugoniot jump 
condition, p{t, 0+) < p < p{t, 0-) and p {S,{t)) < f {p{t, 0±)) for a.e. t € I. Finally, 
by the condition (|4d]l of Definition [U it has to be p {£,{t)) = f {p{t, 0±)) for a.e. t G 
/. D 

The following theorem on existence, uniqueness and stability of entropy weak 
solutions of the constrained Cauchy problem ([T]) is the main result of this paper. 

Theorem 1. Let (F), (W), (PI) hold. Then 

(i) For any initial datum pQ G L°°(]R; [0, i?]), the Cauchy problem ([T|) admits 
a unique entropy weak solution p in the sense of DefinitionUl Moreover, if 
p = p{t, x) is the entropy weak solution corresponding to the initial datum 
po G L°°(M; [0, R]), then for allT > and L > i^ there holds 

\\piT) - P(7')|Ili([-L,L];R) < e IIpo - Po|lLi([-(L+Mr),(L+MT)];R)' (^0) 

where M = Lip(/) and C = 21ji\>{p)\\w\\i^aom_.^y 
(ii) If Po belongs to T>, defined as in ([7]), then the unique entropy weak solution 
of problem ([1]) verifies pit, •) £ V for a.e. t > 0, and it satisfies 

TV {^ (pit))) <Ct = T\ (^ (po)) + 4/(p) + Ct, (11) 

moreover, for a.e. t,s in ]0,T[ we have 

||^(p(t,-)) -^(p(5,-))IIli(R;R) <\t-s\ Lip(*) Ct . (12) 

Proof. The proof consists of three parts, the longest one being postponed to Sec. [3l 

Uniqueness and stability 

Conditions dH of Definition[T]ensure that for all t G [0, T] we can apply the stability 
estimate in Proposition 2.10 in Ref. [1]. More specifically, if p and p are solutions 
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of ([1]) corresponding to the constraints q and g, and the initial conditions po and 
/)o, respectively, then 



\\p{t) - P(0IIl1([-L,L];R) < WPO - ^0|lLi({|x|<L+Aft};M) + ^ / W) - q{s)\ ds. 

J u 

By the explicit expression of the constraints q and g, see ([5]) , we have 



\q{s)-~q{s)\ds= r 
JO 



•p ( / tt;(x) p{s,x) dx] —p [ / ri;(x) p{s,x) dx 



ds, 



and this quantity is bounded by Lip(p) w{0-) f^ \\p{s) - /5(s)||Li([_i,„,o];iR)ds by 
the Lipschitz continuity of p and Holder inequality. Note the inclusions [—L,L] 5 
[— i^, 0] and {|x| < L+MT} I) {\x\ < L+Mi}. We complete the proof by applying 
Gronwall's inequality, see for instance Ref. [29]. 



Existence in V 

The existence problem in the P-framework will be addressed in Sec. O see Propo- 
sition [21 With this result in hand, existence in L°° follows by the density argument 
we develop below. 

Existence in L°° 

Let po be in L°°(R; [0, i?]). By the standard diagonal procedure argument, it is 
enough to prove existence on an arbitrary time interval [0, T] in ]R-|_. Introduce a 
sequence ^q in BV;oc(IK;IK) which converges pointwise a.e. to ^o = '^{po)- Set 
Po = *"H^o)- For any L G N sufficiently large, set p^'^ = p^ X{\x\<l+mt}- 
We have that p^' belongs to V and Pq' converges to po in LJ-^^(M; [0, i?]) as L 
and n go to infinity. Let p"' be the corresponding entropy weak solution of ([1]) 
constructed in Proposition [2l Then, for any L > 1 + i^ 

(A) for all t € [0,T] we have that ||p"'-^(t) — P™''^(*)||Licf_^ n.ig-i goes to zero as 
m and n go to infinity; 

(B) if V > L, then p"'-^ = p"'-^' on [0,r] x [-L,L\. 

Properties (A) and (B) follow by ([10|). Then, by taking L{x) = [|x|J + 1, prop- 
erty (A) ensures that we can introduce the function p(i, x) = lim„_j.+oo p"''^^^\t, x). 
By (B) we also have that p{t,x) = lim„_>._|_oo p"''^ {^ix) for any L' > L{x). 

We prove now that [t >—^ pit)] is an entropy weak solution to ([1]) with initial 
datum pq. For any compact set K C M., take L such that [—L + 1,L — 1] 5 
(K U [-inj,0]). Then p"'-^ converges to p in L^([0,r] x K;[0,R]) and conse- 
quently for a.e. t € [0, T], q^'^{t) = p ( /__; w{x) p'^'^{t,x) dxj converges to 

q{t) = p ( /_j w{x) p{t,x) dx j in L-'-([0,T];M). This is enough to ensure that for 

all test functions (p supported in [0, T] x K, the function p satisfies (|lal) - (|^ and 
that ([5]) holds. In particular the Rankine-Hugoniot condition is satisfied, therefore 
using Lemma [31 we have that f^ (p"' (-,0— )) converges weakly to /(p(-,0— )) in 
Li([0,r];M) and /(p(-,0-)) = /(p(.,0+)). Therefore, also ^ holds true. D 
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3. Wave-front tracking and operator splitting methods 

In this section we construct solutions for initial data in V and we prove that P 
is an invariant domain for the semigroup S. 

Proposition 2. For any initial datum po in T>, there exists a unique entropy weak 
solution of problem ([1]), \t i— )• p{t)] and p{t) belongs to T> for all t > 0. Moreover, 
estimates (fTTj) and (fT2]) are satisfied. 

The solution [t H- p(t)] is the limit (along a subsequence) of a sequence obtained 
by combining the wave-front tracking algorithm and the operator splitting method. 
In the following subsections we describe the construction in full details. 

3.1. Approximation of flux and efficiency functions. Fix /i,n E N suffi- 
ciently large with n^ h. Introduce the mesh 

M^ = f-^2-yimri[o,f{p)]) 

and the set 

P" = PnPC(M;M") , 

where PC (R; M"") is the set of piecewise constant functions defined on M, taking 
values in A4^ and with a finite number of jumps. Approximate the flux / with a 
piecewise linear, continuous flux /" : [0,R] — >■ [0, /(/>)], whose derivative exists in 
[0, R] \ M^ and such that /" coincides with / on M^ , see Fig. [21 left. Clearly, /" 
satisfies condition (F). Consider p'^ {f{M^) np([0, R])) = {|^, . . .,^'^^+2}^ with 
< 1^0 < l^i* < . . . < lm,^+2 ^ -^1 and observe that 

(eti - CT) Lip(p) > P (Iti) - P (^T) = 2- V(P) • (13) 

Approximate p with the function p € PC ([0, R\; /(A4 )) defined as follows: 



/•/" i 



/(M") 




P-P 



/(TW^) 



X" 



RP 



SI ^mi, ^ 



Figure 2. Left: : In thin line / and in thick line the approximation 
/". Right: In thin line p and in thick line the approximation p . 



rrtft-l 



P^iO = Y. Pi X [^h^^h^^ [(0 + Pm, X ^^h^ ^ ^] (0 , (14a) 



CROWD DYNAMICS AND CONSERVATION LAWS WITH NON-LOCAL CONSTRAINTS 







h 
rrih 






h 

rrih- 



R 



where 

^o<^i=ii<---<c, 

PiS.i+i), i = 0,...,mh-l, and p^^ = p(Ch+i) 

see Fig. [21 right. Since h < n, we have that f{M^) C f{M"') 
f{M^') \ {0}, 



p" 



v-p 



L°°([0,/?];IR) 



< r-^f{-p) 



and by ([T3]) and ([T 



inf 

i=0,...,m;j 



P.^ 



t/i 
Si+1 






Si 



> 



2-V(p) , 
Lip(p) 



(14b) 

(14c) 

M]) c 

(15) 

(16) 
(17) 



3.2. The algorithm. Now we can start with the construction of an approximating 
solution \t I— > p'^'^{t)\ to ([T]). As a first step we associate to any fractional time 
interval of the form [^At/j, (£ + l)At/j[, At/j > 0, £ G N, a constrained Cauchy 
problem of the form ([3]) with constant constraint. Then the wave-front tracking 



n,h 



algorithm gives us the corresponding exact solution [t !->■ p^]^^{t)\. Finally, p"' is 

obtained by gluing together /o"^]^, ^ € N. The existence of a limit for p^'^ as n and 
h go to infinity is ensured by the choice 

^^^ ^ 2^+iu;(0-)Lip(p) ' ^^^^ 

which will be motivated in the proof of Lemma [1] by a sort of CFL condition. 
Roughly speaking, this condition is needed to bound the possible jump in the 
value of the constraint due to the update at each fractional time (£ + l)Atft, £ G N. 
Approximate po with a piecewise constant function pg : M — )• [0, R\ that coin- 
cides with pq on M.'^ and such that II/Oq IIli(r-r) — II/'oIIlicr-r) ^"^^ TV (^ (Po)) ^ 
TV(^(/9o))- Clearly, p^ belongs to P". First consider the approximating con- 



strained Cauchy problem 

dtp + d^r{p) = {) 

r(p(t,o±))</(ss) 

p(0,x)=p^(x) 
where 



(t,x)G]0,Aift[x 
iG]0,Aife] 

xGM , 



"0 



w{x) Po(x) dx . 



n,h/ 



The unique exact solution [t >-^ p^' (t)] for the above problem is obtained by piecing 
together the solutions to the Riemann problems at points where p^ is discontinuous 
or where interactions take place, namely where two or more waves intersect, or 
one or more waves reach x = 0. For the definition of solution of the Riemann 
problem with a piecewise linear, continuous flux away from the constraint, we 
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refer to Sec. 6.1 in Ref. [2] or to Sec. 5.2 in Ref. [29]. The definition of solution to 
the constrained Riemann problem along x = follows by the obvious adaptation 
of Definition 2.2 in Ref. [8J to the case with a piecewise linear continuous flux, 
see also Sec. 6.3 in Ref. [29J. The results of Theorem 3.4 in Ref. can be easily 
generalized to the case with piecewise linear continuous flux and, therefore, we can 
define 

/o"'''(t, x) = p'l'^it, x) for (t, x) G ]0, Ath] x M . 

We can assume that no interaction occurs at time t = Ath, see assumption H2 
below. Then the approximate solution is prolonged beyond t = Ath by taking 

p'^'^'(t, x) = pl'^{t - Ath,x) for (t, x) G ]Ath, 2Ath] x R , 

where [t i— ?> p^' (t)] is the exact solution of the constrained Cauchy problem 
dtp + d,rip) = {t,x)e]0,Ath[xR 

r(/,(t,0±))</(H^'") te]0,Ath] 



with 



p{0,x) = p"' {Ath,x) X G 



^"' = / w{x) p"' {Ath,x) dx . 



We repeat this procedure at each fractional step and, once we get [t i-^ p"' (t)], 
we construct [t i-^ Pilii^)] by solving a constrained Cauchy problem of the form 

dtp + d^r{p) = {t,x)e]0,Ath[xR (19a) 

r {p{t, 0±)) < / (h"''^) t G ]0, Ath] (19b) 



piO,x) = p';'\Ath,x) xGM, (19c) 



where 

^n,h I / \ n,h 



I w{x) pl^{Ath,x)dx . (19d) 



We stress that the solution to (|19p is unique and that the efficiency at the exit 
may change at each time t G At/^N and only there. 

To simplify the wave-front tracking algorithm, see Remark 7.1 in Ref. j2], it is 
standard to remark that, without loss of generality, one can assume that: 

HI At any interaction either exactly two waves interact, or a single wave reaches 

the constraint x = 0. 
H2 No interaction occurs at time t G At/jN. 

In this way we construct 
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and an approximate solution of the Cauchy problem ([T|) 



p-'\t, x) = Y, p,ti(^ - ^^t'^^ ^) XuAt,, (i + l)At J (*) ' (21) 

where [t i->- p^!^_i{t)] is the unique solution to (fT9|) . 

Roughly speaking, the present procedure consists in the application of two op- 
erators, and S, at each fractional step ]iAth,{i + l)Ath], i & N. The first 
"' = Q[p^' (At/j)], while the second operator gives the solution 
S[p^' (At/j),H"' ] of the constrained Cauchy problem of the form (fT9|) . 

with [x I— )• yo"' (At/j,x)] as initial datum and with p''(H"' ) as constraint. 
More rigorously, for any pg G 2?" and t G M-|_, define recursively 



operator gives 

n,h 
Pi+l 



if t G [0,Ath], and, if t G ](^ + l)At/j, i 



+ 2)At/,],^GN, then 



-in,/i 



F"''^K](t) = 5 F"''^[pS]((^ + l)At,),e 



nn,/l 



:pS]((£+i)At,) (t) 



3.3. ^ priori estimates. In this section we prove that p^'^{t) = F^'^[pQ][t) is in 
D" on any bounded time interval [0,T], T > 0, and we estimate TV (^ {p^'^{t))) 
uniformly in n, h and t. To this aim, we introduce the following Temple functional 



with 



T^''' (t) 



7'(P,H) 



Tv('f (p"''^(t))) +7^^ (p"''^(t),H"''^(t)) +r5^(t) , (22) 





A[f{p)-p\'B)\ 
" T 



if p(O-) > p > p(0+) and 



r(Ko±)) 

otherwise, 



p" 



m) 



r!i.{t) = 5-2-^ fip) 



At. 



t 



Ah 



where [-J : M — )• Z denotes the floor function. Recall that the Temple functional 
adopted in Ref. [8j involves the total variation of the approximating constraint, 
which is given a priori. In our construction, at each fractional time interval we 
are dealing with a different approximating problem ()19p and we need to know the 
solution at the previous step in order to fix the value of the constraint in (]19bp . 
Therefore, the constraint p (H"' (t)), t G M+, and its total variation are not given 
a priori. Nevertheless, due to the choice of Ath, we are able to bound the possible 
jump of p (H"' (t)) at each time step, as we will see in Lemma [H and estimate a 
priori the total variation of the efficiency. From this point of view, the functional 
T^' is the natural generalization of that one used in Ref. [8j. In fact, the two 
functionals have in common the first two terms, namely 



Q"''^(t)=Tv('f (p"''^(t))) 



^h fn.h 



{t),E-''{t)) 



while r^ (t) is introduced to control the total variation of p (E"''^( 
interval [t,T]. 



(23) 
■)) in the time 
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Lemma 1. For any £ G N, the jump in the efficiency at time t = {i + l)Ath, 
namely 



hf'^n,h\ „hf'=n,h 



p (H^' -^) — p (H^' ) , is either zero or 2 f{p). 



Proof Fix ^ G N. If 






?i+l Si 



then [^ iH^ V^^K)\ has 



< inf 

i=0,...,r?i;j 

at most one jump for ^ between H"^-,^ and H"' and ([T6|) aUows us to conclude. 
Because of (fT7|l . we just need to show 






Lip(p) < 2-V(p) . 



By Proposition[Tl p^' ^ is a weak solution of the problem (|19ap . (|19cp with p^' (At/j) 
as initial condition. Then, for any (p in Cl.{M?;'. 



we have 



n,h 



n,h 



n,h/ 



Pt+i dt(l} + f {p^i+i) d^(j) dx dt+ / p"' iAth,x) (/'(0,x) dx = . (24) 



Let {r]u)i, be a standard family of mollifiers and define w^, = w * rj^. Let de be as 
in Q. Take < ti < t2 < At/j and consider the test function 



(t>{t,x) 



t-t-i 



5e{z) dz 



t~t2+e 



^eiz) dz 



X+£ 



Wu{x) 



Observe that (/>(0, •) = and that letting e go to zero we get 



We pass to the limit in the Eq. (j24p letting e go to zero and we obtain 

dx 



Wu {x) p"|i (ti , x) - p"|^ (t2 , x) 



t2 



Wu{Q-) /(/5"+i(t,0-)) -Wu{-iw+) f [p'l:!^^{t,-\n, + ) 



dt 



Then, as v goes to infinity we get 

w{x) p^^'^^ (ti , x) - p1'_!^^ (t2, x) 

) 

By (fT9c]l and (fT9dl> we have that 



dx 



< (t2 - ti)/(/5M0-) . (25) 






w[x 



/5"+i(Ath,x) 



n,h 



n,h 



{Ath,x) 



dx 



n,h 



w{x) p"_;^(At/j,x) -p"|^(0,x 



dx 



< Ath f{p) u)(0- 



Therefore by (jlSp the proof is complete. 



D 



n,/i • 



We are ready to show that T^' is a Temple functional 
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Proposition 3. Let /i,n G N and p^ E D". On [0,T], the map [t ^ T^'^ (t)] is 
non-increasing and it decreases by at least 2^"f{p) each time the number of waves 
increases. 



The proof is deferred to Sec. 16.11 

In the next corollary we rely on Proposition [3] to prove a uniform estimate on 
TV(^(/j"''*(t))). 

Corollary 1. There exists a constant C > 0, that does not depend on n or h, such 
that for all t > 

TV (^ (p"'^t)) ) < TV (^ (po)) + 4/(p) + Ct. (26) 

Proof. We consider the functional Q"''^ = T^' — Fj. introduced in (j23p . Proceeding 
as in the proof of Proposition [3l we can show that Q"'^ may increase only at 
t € AtfiN. However, since T^' is strictly decreasing at t G At/jN, we have that for 
all ieN, 

Q"''*(Mt/,+) - Q'^'^iAth-) < T!^{iAth+) - rUeAth-) = 5 • 2-^f{p) . 



Therefore, by UE 



TV U- p"''^(t) < Q"'^(t) < Q"''^(0) + 5 • 2-V(/o) 



t 



< TV (^ (po)) + 4/(p) + 10 w{0-) Lip(p) f{p) t , 
and the estimate 1^ holds with C = 10 w{0-) Lip(p) f{p). D 

By the results proved in Ref. [8], the assumption H2 and the corollary above, 
we have that both p'^'^{t) and rP'^{t) are well defined for any t E [0, T] and that 
pn,h ijgiongs to C'^(M+;P'^). In particular, [t i— >• p^'^^t)] is piecewise constant 
with discontinuities along finitely many polygonal lines with bounded speed of 
propagation, that do not intersect each other at any time t E At/jN. By the 
construction of p^'^ and its continuity with respect to time it is not difficult to 
show 

Proposition 4. The map \t i— )■ p"' (t)] given by ()2ip is an entropy weak solution 
in the sense of DefinitionU\ (with f^,p^ replacing f,p) to the problem 

dtp + 9^/"(p) = (t, x) E M+ X M (27a) 

r {p{t, 0±)) < / (H"'^(t)) t E M+ (27b) 

p{0,x)= p'Six) xeR, (27c) 

where [t \-^ H"''^(t)] is given by ([2n|) . 

The proof is deferred to Sec. 16. 2[ 

Proposition 5. There exists a subsequence of p'^' converging a.e. on M+ x R to 
a limit p E L°° (M-|_ x M; [0,i?]). In addition, p satisfies estimates (jlip and (|12p . 
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Proof. By the standard diagonal procedure argument, it is enough to prove con- 
vergence on an arbitrary time interval [0,T], T > 0. The sequence ^ (p"' ) is 
uniformly bounded in L°° ([0,r] x M;M) n L~ ([0, r];BV(M;R)) by Corollary [H 
In order to get compactness in LJ-^^, see Theorem 2.4 in Ref. [2], we still need to 
show that [t 1-^ ^ (^p"''^{t, •))] is Lipschitz with respect to the L-'--norm. In analogy 
to (0, we define 



^^ip) = signip-p) irip)-r(p)) 



rir) 



dr . 



(28) 



We observe that \I' coincides with ^" on M^ and, as a consequence, Lip(^"') < 
Lip(^). Since p"-''^ takes values in 7W", we have *" (p"-'') = ^ (/>"•''). Hence, 
TV (^" (/o'''''(t))) = TV (* (/9"'''(t))) and by Corollary [U we have 



n I ^n,h 



5..^" p^ 



L-^{[0,T];Mbi 



<CT = TV('f (po))+4/(/?) + C7r 



uniformly in n and h. Above, A^fc(M;M) denotes the space of bounded Radon 
measures. Let 5" = /" o (^")~^ and remark that by (j28p 



^n(^) ^ r°(^")-HV') ^ roi^n-'w ^ {_^^ i| _ 



^"o (^")-l(V') 



Hence dtp""^' is bounded in L~ {[0,T];M 
rem 4 in Ref. |2I] . we have 



/"o(1'")-l(^) 

because, by Eq. ([?7|) and Theo- 



Stp"-'^ 



L°°([0,r];A1i,(K;I 



< ll5"llL-{h/(p),/(p)];I 



a^^" ' ""''' 



L°°([0,T];A1i, 



As the functions \I'" are uniformly Lipschitz, also the distributions /u"'^ 
dt^'"{p"''^) are uniformly bounded measures in L°^ ([0, T];7W;,(M;M)) with 



,n,h 



fJ' 



L°°([0,T];A^5(IR;]f 



< Lip(^") Ct ■ 



Now, let {t]u)u be a standard family of mollifiers in C^(M^;! 

q,n ^pn,h^ ^ ^^ ^j^j ^n, _ ^n,h ^ ^^_ Then 



and define F^ 



n,h 



Ml 



< 



/^ 



,,/i 



L°°([0,T];A46 



L°°([0,T];Li 

Due to the regularity of Fy' , for any 5 > and for any 0<t<t + (5<T 

Li(IR;IR) J 

<6 



\n,hi 



F;^'''{t + 5,-)-F^'\t 



l^T{s,x)ds 



n,h 
Pu 



L°°([0,T];Li(R;I 



dx 

< 6 Lip(*") Ct , 
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and as v go to zero we deduce the uniform Lipschitz continuity in time of ^ (p"' ) = 

^{p-'\t + 5, •)) - ^{p'^'Ht, •)) , ,,„ „, < 5 Lip(*) Ct . 

In this way we prove the existence of a subsequence of ^ (p"''^) = 'I'" (p"''^) that 
converges in Lj^^{[0, T] x M; M) to a function V in L°° ([0, T]; BV(M; [- f (p) , f (p)])) 
which satisfies 

mt + 5,-)- i'it, OIIlicr^ < 5 Lip(*) Ct ■ (29) 

For simpHcity we stih denote the subsequence ^ (p"' ) . Since ^ is invertible 
and '^~^ is continuous, also p"''* converges in Ljjj^([0,T] x M;M) to a function p = 
*-i(V') in L°° ([0, T] X M; [0, i?]). In particular, by (I2SD and (I22D the estimates (HI]) 
and (USD hold true. 

Lemma 2. For any T > 

r-T 



D 



lim 



H"''^(t)- / u;(x) p{t,x) dx 



dt = . 







£=0 



£At. 



Proof. Let T > and define ^^ = [T/Ath] ■ Then by i^ and 
S"''^(t)- / t/;(x) p"'''(t,x) dxdi< 

R_ 

Ef'^{t - tAth) - / w{x) p"l^i(i - Mt/i, x) dx 
+ /' H;^"(t - 4At,) - / «;(x) p;^" (t - 4At;„ x) dx 

^(a^) Pm (0' 2;) - p";'\ (t, x) 

/(p) ^(0-) 



dt 



dt 



T rAth 





E 

£=0 

+ 



dx 



dt 



< 



T-^Ath 





V^ 1 

T rAth 



dx 



dt 



t dt + 

^0 



E 

Atl 4 + (T - ^At,,) 



T-^Ath 



t dt 



/(p) u;(0- 



Therefore, since if^Ath converges to T as h goes to infinity and p"''^ converges to 
p in LJ-^p(]R-|_ X M;R) as n and /i go to infinity, the proof is complete. D 

Since p"''^ converges to p in Lj-^^, Proposition U] and Lemma [2] imply that 
[t I— >• p(t)] satisfies the conditions ([la|) -([Ic l) and ([5]) of Definition [1] with respect 
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to the problem ([1]). Moreover, p satisfies the condition (|4dp of Definition [1] by 
Lemma [3l and it satisfies estimate (fTT]l and (fT2]l by Proposition [5l Finally, ob- 
serve that p G C'' (M+;LJjj^(M; [0, i?])) because of entropy inequalities (fla|) - ()ic|) . 
see Ref. [3J. As already observed ^ {p) G C° (E+;BV(M;M)) thus p{t) G V for ah 
t. 

Uniqueness of the entropy weak solutions to the Cauchy problem ([1]) in the case 
p G Lip([0, i?];M), po (zT) follows directly from uniqueness in the L°°-framework, 
see the first part of the proof of Theorem [TJ 

4. The constrained Riemann problem 
In this section we study constrained Riemann problems of the form 

dtp + d^fip) = (t, x) G M+ X M (30a) 

f{p{t,0±))<p( w{x) p{t,x) dxj teR+ (30b) 

P{0,X) = { ^^ -r"^ Jn XeW (30c) 

[ PR 11 X >0 ^ ' 

with PL, PR € [0,^]- Along with (F) and (W), we assume that: 

(P2) p belongs to PC ([0, R] ; ]0, /(/>)]) and is a non-increasing map. 

The assumption (P2) is introduced in place of (PI) to allow an explicit con- 
struction of solutions to (f30|) . However, the regularity of p required by (P2) is 
not enough to apply the results of Theorem [T] In fact, the uniqueness of entropy 
weak solutions as well as the stability estimate (|10p do not hold in the present 
framework, as we will see in Example [2J 

Aiming for a general construction of the solutions to (|30p . we allow p to be a 
multi-valued piecewise constant function, namely, see Fig. [TJ right: 

• there exist S,i,... ,S.n G ]0,i?[ and po,- ■ ■ ,Pn G ]0, fip)], with S^i < ^j+i and 
Pi > pi+i, such that p{0) = po, p{R) = Pn, p Xit. t. r = Pi for i = 0, . . . , n, 

p{Ci) = \Pi,Pi~i] for i = 1, . . . , n, being Co = and C„+i = R. 

Let a{pL,PR) = U{pl) - f{pR))/{pL - Pr) be the speed of propagation of a shock 
between pi and pR, while \{p) = f'{p) is the characteristic speed. Introduce the 
maps p,p : [0, /(/?)] — )• [0, i?] implicitly defined by 

/ {p{p)) =P = f {p{p)) and p{p) < p< p{p) . 

Let Pi = p{pi) and pi = p{pi). Denote by TZ the classical Riemann solver Refs. [55], 
[25] . This means that the map [(t, x) >->■ TZ[pL-,PR\{x/t)\ is the unique entropy weak 
solution for the unconstrained problem (|30ap . (|30cp . see Refs. [2], [2S], [3D] for its 
construction. As we will see in Proposition [UJ the classical solutions given by TZ 
may not satisfy the constraint ()30bp . For this reason we consider also nonclassical 
solutions, namely solutions that do not satisfy the Lax entropy inequalities, see 
Ref. [24J as a general reference. In general, entropy weak solutions to ()30p are not 
self-similar nor unique, as we will show in the two following examples. 
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Example 1. Let < pL < pR < R be such that fipi) > Kpr)- If Ci < Pl < ^i+i, 
Pi+i < fipn) < Pi and j > i is such that f (pj) = p (pj) and f (pk) > p (pk) for all 
k € {i, ■ ■ ■ , J — 1}, see Fig. left, then the entropy weak solution to the correspond- 
ing Riemann problem (J3U|) is not self-similar, see Fig. 0, right. More in detail, for 
sufficiently small times, the solution corresponds to the classical one and is given by 
a shock with speed a{pL,pR) < 0. The corresponding map [t i— t- ^{t)] is increasing 
and by hypothesis, there exists a time ti < —iw/cr{pL, Pr) such that ^(ii) = Ci+i- 
Then, the efficiency of the exit falls to Pi+i and the solution given by the classical 
Riemann solver TZ no longer satisfies the constraint condition ()30bp . As a result, at 
time ti the solution performs a nonclassical discontinuity at the constraint location 
and two further classical shocks appear, one with speed a{pji, pi+i) < and one 
with speed a{pi^i, pr) > 0. The final solution can then be constructed by taking 
into account the interactions between the shocks on each side of the constraint and 
the appearance of new shocks each time [t H' ^(t)] crosses ^k, k € {i + 1, ■ ■ ■ ,j — l}. 





PR Pj P 



Figure 3. Construction of a non self-similar entropy weak solution 
as in Example [TJ On the left, the thick line corresponds to the 
efficiency of the exit p\]p^^p[. 



As we have seen, the lack of self-similarity is related to the jumps of [t i— )• 
p(^(t))]. Nevertheless, in the proof of Proposition [6] we show that any entropy 
weak solution of ()30p is self-similar for sufficiently small times. Therefore, it makes 
sense to introduce nonclassical local Riemann solvers, see Definition [2j Then, the 
availability of a local Riemann solver allows us to construct a global solution to 
the Riemann problem ()30p by a wave-front tracking algorithm in which the jumps 
in the map [t \-^ p{^it))] are interpreted as interactions. 

The next example shows that the entropy weak solutions to the constrained 
Riemann problem ()30p are not necessarily unique. 



Example 2. Consider the constrained Riemann problem (j30p with pi = S,i+i € 
]p,R[ and pR = p. Assume that f {pi+i) = Pi+i < /(Ci+i) < Pi = f {pi) < f{p), 
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Pi+l Pi PiPLPi+lP 



Figure 4. With reference to Example [21 the flux configuration 
and three different solutions pi, p2 and p^ to the same Riemann 
problem are represented from left to right. Here pi = p{f{pL))- 



see Fig. [7| left, then 



piix/t) 
P2{x/t) 



T^[Ci+i,Pi+i]ix/t) ifx<0 
T^[Pi+i,p]{x/i) ifx>0, 

Ci+i ifx<0 

n[p{f{^i+i)),p]{x/t) ifx>0 

TZ[^i+uPi]{x/t) tfx<0 
TZ[pi,p]{x/t) ifx>0, 



are self-similar entropy weak solutions of problem (|3Up with the same datum, see 
Fig.\^ Clearly, the above solutions are distinct if Pi+i 7^ /(Ci+i) / Pi; otherwise 
two of them may coincide. Additionally, for an arbitrarily chosen t > 0, the 
functions 



Pt,i{t,x) 



PtAt,x) 



P2 (x/t) ifO<t<t 

pi {x/{t — t)) if t > t and a; < 

Pi+i ifi_<t<ii andO <x< a {pi+i, p{f{Ci+i))) {t - t) 

p ift<t<tiandx>(j{p{f{ii+i)),p)t 

Pi+i if t > ii and < x < a (/Oj+i, p) {t — ii) + xi 

p if t > ii and x > a {pi+i,p) [t — ii) + xi 

p{f{ii+i)) otherwise, 

P2 {x/t) ifO<t<i 

p3 {x/(t — t)) if t > i and x < 

Pi ifi_<t<ij, andO <x< a {pi, p{f{(,i+i))) {t - t) 

p ift<t^<izandx>(j{p{f{ii+i))^p)t 

Pi ift> is and < x < a (pi, p) {t - i^) + X3 

p iftyi^ and x >a {pi, p) {t - is) + £3 

Hf{^i+i)) otherwise, 
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where 



h 



t o-(/5i+i,p(/(^i+i))) 



h 



a ipi+i,p ifiCi+i))) - o- (/5 (/(Ci+i)) , p) 
i a{pi,p{f{ii+i))) 



xi =ii cr(/5(/(^i+i)),p) 



^3 = ^3 cr(/5(/(^i+i)),p) 



£^ {pi,P (/te+i))) - 0- (P (/(Ci+l)) , P) ' 
are also entropy weak solutions, see Fig. O Therefore, because of the arbitrariness 
of t, we can build infinitely many different solutions which are not self-similar 
on any open time interval. However, remark that the asymptotic profile of pn 
coincides with the asymptotic profile of pi . 

One may guess that the lack of uniqueness is due to the fact that p is a multi- 
valued function. But we observe that if we pick up /(^j+i) as value for p{S,i+i) , then 
still all the above solutions are admissible. Moreover, if we have p{S,i+i) ^ /(^j+i), 
then both pi and p^ are admissible, but not p2- Thus there is more than one entropy 
weak solution even if p is a single valued function. 





Figure 5. With reference to Example [21 the solutions p^i and p^^. 
Above PL = p{f{pL))- 

Introduce the subset of [0, R] 

C = Upl,Pr) G [0,Rf ■■ (pL,pfi) satisfies condition (C)| , 

where we say that {pl, Pr) satisfies condition (C) if it satisfies one of the following 
conditions: 

(CI): PL < PR, fipR) < f{pL) and f{pR) < p{pL+); 

(C2): PL < PR, fipL) < f{pR) and f{pL) < p{pl+); 

(C3): PR<PL<P and /(p^) < P {pL+); 

(C4): pR<p< PL and f{p) =p{pL+); 

(C5): p < PR< PL, f{pR) <P{PL-) and /(pl) <p{pL+)- 
In Proposition [6] we will prove that a constrained Riemann problem admits as 
unique entropy weak solution the classical one, at least for small times, if and only 
if its initial datum satisfies condition (C). 
Analogously, introduce the subset of [0, R] 

J\f = Upl,Pr) e [0,Rf ■■ (/?L,/5i?) satisfies condition (N)| , 

where we say that (pl, Pr) satisfies condition (N) if it satisfies one of the following 
conditions: 
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(Nl): PL < PR and /(/jl) > /(p^) > PipL+); 

(N2): PL < PR, f{pL) < f{pR) and /(pl) > pipL-); 

(N3): PR<PL<P and /{pl) > P (pL-)] 

(N4a): pR < p < pl, f{p) y^p{pL-) and /(pi) < p{pL+); 

(N4b): pR<p< PL, f{p) 7^p{pL-) and /(p^) > p(pl-); 

(N5a): p < pr < pl, /{pr) > P (pl-) and f{pL) < p {pl+); 

(N5b): p < PR< PL and /(p^) > p [pL-)- 
In Proposition[6]we will prove that, at least for small times, a constrained Riemann 
problem has a unique entropy weak solution which is nonclassical if and only if its 
initial datum satisfies condition (N). 

Observe that if the constraint function p is constant in a neighborhood of the 
state PL, then p{pL—) = p{pl+) and this simplifies the above conditions. Also 
a right or left continuity assumption on p would simplify the above conditions. 
However, we keep p as a multi-valued function to take into account all the possible 
solutions of (pOj) . 

In the next proposition, we show that uniqueness holds if and only if the initial 
data are in C U M. 



Proposition 6. Consider the constrained Riemann problem ()30p . 

• // {pl,Pr) € C, then the map [{t,x) i-^ TZ[pL, PR]{x/t)] is the unique entropy 
weak solution at least for t > sufficiently small. 

• If {pl,Pr) S J^ , then there exists a unique p E [p{pL+),p{pL — )\ such that the 
map 

' { n[pL,p{p)]{x/t) ifx<o' 
\n[p{p),pR]{x/t) z/x>o 

is the unique entropy weak solution at least for t > sufficiently small. 

• // (pl, Pr) G [0,-R]^ \ (C U Af), then the corresponding constrained Riemann 
solver ()30p admits more than one entropy weak solution. 



Proof. We stress that any nonclassical entropy weak solution in the sense of Def- 
inition [T] is also a classical entropy weak solution in the Kruzkov sense in the 
half-planes M+ x M_ and M+ x M+. Therefore, at least for t > sufficiently small, 
by Proposition [T] and assumption (P2) any nonclassical entropy weak solution 
of PU|) must have the form, see Fig. EJ 

[ n[p{p),pR\{x/t) ifx>o. 

Observe that ()31ap is uniquely identified once we know p which, by ([8]), satisfies 

P = f{p{p)) = f{pip)) ■ (31b) 

We recall that (|31bp means in particular that the Rankine-Hugoniot jump con- 
dition is satisfied at x = even when the solution to the Riemann problem is 
nonclassical. As a consequence of (j31bp . of assumption (P2) and of the continuity 
of [t I—)- C(t)], we have that 

pe [p{pL+),p{pL-)]- (31c) 
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Figure 6. The four possible configurations of nonclassical entropy 
weak solutions of the form (j31|) . 



This implies that if p{pL+) = p{pl—)-, then p(^) is constant in a neighborhood 
of PL and, since the solution is in C° (M+;Lj-^p(IR; [0,i?])), uniqueness is ensured 
by the results in Ref. [8]. However, the continuity of p at pL is not a necessary 
condition for uniqueness. In Sec. 16.31 we prove that: 

{pl,Pr) G C: In this case, the corresponding classical solution satisfies (pO]) 
for all t > sufficiently small and it is not possible to construct a different 
solution. 
{pl,Pr) G a/": In this case, the corresponding classical solution does not sat- 
isfy (jSObp , and there exists a unique nonclassical solution that satisfies (j30|) . 

Now we list the "pathological" cases, where we have more than one admissible 
solution. We stress once again that a necessary condition for non-uniqueness is 
p{pl—) 7^ p{PL+) and p{pL—) > fipi) > p{pl+)- It is important to stress that in 
general the solutions to the constrained Riemann problem (|3Up are not self-similar, 
see Example [H All the cases listed below describe self-similar solutions because 
we let the solutions evolve only on a small interval of time. 

(CN2): If PL < PR, f{pL) < Upr) and p{pL+) < /(pl) < p{pl-), then 
the classical solution [{t,x) i-^ 7l[pL, PR\{x/t)], which consists of a shock 
with non negative speed, as well as the nonclassical solution pip , with 
P = p{PL+), are distinct solutions of (pO]l . 

(CN3): If PR < PL ^ P and p{pL+) < Upl) < p{pl—)-, then the classical 
solution [{t,x) i-7> IZIpl, PR\{x/t)], which consists of a possible null rarefac- 
tion on the right of the constraint, as well as the nonclassical solution ()3ip . 
with p = p{pL+), are distinct solutions of (pO]l . 

(NNN4): If PR < p < PL, P{PL-) 7^ P{PL+) and p{pL+) < /(pl) < 
p{Pl — ), then the nonclassical solutions of the form (f3T]l corresponding 
to p G {p{PL+), f{pL),p{pL—)} satisfy PU|) . see Example [2J Observe that 
such solutions are distinct as far as they correspond to distinct constraint 
levels p, and that in any case there exist at least two distinct nonclassical 
solutions. 

(CNN5): If p< PR < PL, Kpr) <p{pl-),p{pl-)^p{pl+) andp(pL+) < 
f{pL), then the classical solution [{t,x) i-^- TZ[pL, Pr\{x /t)], which consists 
of a possible null rarefaction on the left of the constraint, as well as the non- 
classical solutions of the form (f3T]) corresponding to p G {p{PL+), fiPL)} 
satisfy (|30p . Observe that the two nonclassical solutions are distinct as far 
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as they correspond to distinct constraint levels p, and that in any case there 
exist at least two distinct solutions, one classical and one nonclassical. 
(NNN5): Up<pR< pL, /(pr) > p{pL-) > Kpl) > p{pl+) andp(pL-) ^ 
p{pL+), then the nonclassical solutions of the form (|31|) corresponding to 
p S {pipL+), fiPL),p{pL—)} satisfy ([5U|) . Observe that such solutions are 
distinct as far as they correspond to distinct constraint levels p, and that 
in any case there exist at least two distinct nonclassical solutions. 

This concludes the proof. D 

As the local solutions of the Riemann problem are not unique in general, we 
are naturally led to question the existence of suitable selection criteria. All the 
solutions we introduce are solutions in the Kruzkov sense in the open half-planes 
M_i_ X M_|_ and M+ x M_, so they satisfy the basic requirement of entropy dissipation. 
However, coming back to the real situation which our model aims to describe, 
namely the evacuation of a narrow corridor, we argue that the most desirable 
solution is obviously the one corresponding to the highest admissible values of the 
flux at the exit. In analogy to the discussion in Ref. [19j we interpret all other 
possible solutions as consequences of an irrational behavior, which in literature is 
often described as panic. It is also important to remark that since non-uniqueness 
is possible only when p{pL-) 7^ p{PL+) and p{pL-) > fipi) > p{pl+), non- 
uniqueness concerns at most a finite number of left states. 

From now on we restrict ourselves to the case in which p{£,i) can only take the 
values Pi and Pi+i and not the intermediate values, because the extremal behaviors 
are the most relevant in view of the applications. 

Definition 2. Two Riemann solvers IZ'^ and TZ^ for (j30p are defined as follows 
for t > sufficiently small and a; € M; 

(C): If{pL,PR)eC then 

n'^[pL,PR]it,x) =nP[pL,PR]it,x) =n[pL,PR]ix/t). 

(N): If{pL,PR) eAf then 

' 7^[pL,/5(p)](x/t) i/x<0 



where p = p{pL—) if {pl,Pr) satisfies (N^a) or (N5a), otherwise p = 
p{pL+)- 
(CN2), (CN3), (CNN5): // {pl,Pr) satisfies one of these sets of condi- 
tions then 

T^''[PL,PR]{t,x) =n[pL,PR]{x/t) , 



nP[pL,PR]it,x) 



n[pLrp{p{pL+))]{x/t) ifx<0 
n[p{p{pL + )),PR]{x/t) lfx>0. 



(NNN4), (NNN5): If(pL,pR) satisfies one of these sets of conditions then 
TZ'^Ipl, PR\{t,x) takes the form ([3T]l with p = p{pL—) and TZ^Ipl, pR]{t,x) 
takes the form (f3T]l with p = p{pL+)- 
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In the next proposition we collect the main properties of the Riemann solvers 
TZ'^ and TV . In particular (R6) means that the Riemann solver Tf^ is the one 
which allows for the fastest evacuation, while TW is associated to the slowest one. 

Proposition 7. Let {pl,Pr) G [0,R] ■ Then, for -k = q, p: 

(Rl) [{t,x) 1-^ TZ*[pL-, PR\{t.,x)\ is a weak solution to ()30ap . ()30cp . 
(R2) TZ*[pL,pfi] satisfies the constraint ()30bp in the sense that 

f{n*[pL,PR]{t,0±))<p( [ w{x)n*[pL,PR]{t,x) dx 

\Jr- 

(R3) n*[pL,PRm e BV(M; [0,R]). 

(R4) The map 7^* : [0,R]'^ -^ Lj-^^(M+ x R;R) is continuous in C U Af but not 

in all [0,i?]2. 
(R5) TZ* is consistent, see Refs. [8], [12]. 
(R6) TZ'^[pl,pr\ maximizes the flux at the exit, in the sense that if £ is the set 

of all entropy weak solutions of the Riemann problem, (|30p . we have 

max{/(p(t,0±))} = /(7^'?[/,L,/0K](0±)) . 

Analogously, TU'lpi, Pr\ minimizes the flux at the exit, in the sense that 
min {/(p(t, 0±))} = / (7^P[/,L, pr]{Q±)) . 

pGt 



The proof of Proposition [7] is deferred to Sec. 16.41 

It is important to observe that even if p{^i) can only take the two values pi 
and Pi+i, this is not enough to rule out the existence of infinitely many different 
solutions as the ones described in Example [21 in the case pi > f{Ci) = Pi+i- 
However, the Riemann solver TZ* spontaneously selects one of them because it 
sticks to the constant level of constraint prescribed by Definition [2] until a nonlocal 
interaction takes place. 

Remark 1. Although the Riemann solvers TZ* are not Lif^^-continuous, an exis- 
tence result for the Cauchy problem ([T|) can be obtained from a wave-front tracking 
algorithm based on TZ* , see for instance Ref. |16j. Such approach using TZ* does 
not require the operator splitting method. However, the non-local nature of the ap- 
proximating problems prevents us from a direct application of the Riemann solvers 
IZ* . In fact, even in a arbitrary small neighborhood of x = 0, to prolong the 
approximating solution p" beyond a time t = i > it is not sufficient to know 
the traces p"{i,0—), p'^{t,0+), but also the value /_. w{x) p'^{i,x) dx is needed. 
Roughly speaking, because of the non-local character of the constraint one cannot 
merely juxtapose the solution to the Riemann problem associated to the values of 
the traces at x = with the solution to the Riemann problems away from the con- 
straint. Finally, also jumps in [t i— )• p{^{t))] have to be considered as (nonlocal) 
interactions. Therefore, the approach using TZ* is considerably heavier and more 
technical than the one we presented in Sec. \3.iA and we do not pursue this line in 
this paper. 
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(a) The solution in tiie (a;, i)-plane. 




(b) Profiles of [x i-> p{t, x)] at 
times t = 0, £b/2, {tp + £g)/2, 
(tG + tj/)/2, (tQ+ifl)/2. 



Figure 7. The solution described in Sec. [5] and corresponding to 
the choice (1321) 



5. Numerical examples 

In this section we apply the model ([T]) to simulate the evacuation of a corridor 
through an exit placed in x = 0. The simulation is obtained by explicit analysis 
of the wave front interactions, with computer-assisted computation of front slopes 
and interaction times presented on Fig. [71 

Assume that the pedestrians are initially uniformly distributed in x E [xa,2;_b] 
with maximal density, namely po = R X[xa,xb]- ^^ ™ -^^S- El (a), we choose the 
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efficiency of the exit, p, of the form 

( po if < e < 6 

p(0 = { Pi if 6 < C < ^2 

( P2 iH2<C<R, 

and such that the solution to each Riemann problem is unique and TZ^ = IZ'^. 

Then we can start with the construction of the solution. From B = (xb,0) 
starts the rarefaction TZb that takes the values TZB{t,x) implicitly given by 

X{R) < A(7^B(t,x)) = ^^ < A(0) . 

The first pedestrian reaches the exit at time tc = — xb/A(0). In L = (xaj^l); 
with ti = (xA — xb)/X{R), the stationary shock Ca originated from A = {xa,0) 
starts to interact with the rarefaction TZb- As a result, from L starts a shock Cl 
given by 

Cl- x{t) = a{0,TZBit,x{t))) , x{tL) = xa- 

At time t = to the maximal efficiency of the exit pq is reached, / (7?.B(t£), 0)) = pO) 
and a queue appears behind it. The tail of the queue is represented by the backward 
shock Cd given by 

Cd: x{t) = a{nB{t,x{t)),p{po)) , x(tD) = 0- 

Let tE be the value of t solving the equation 

E: / w{x) TZB{t,x) dx + p{pq) / w{x) dx = ^i, 

t > to , Coit) > -itu . 

The data can be chosen in such a way that CiitE) < —iw, then at time t = Ie 
the efficiency of the exit falls to pi and a further shock Ce with constant speed 
a {p{po) -, p{pi)) < appears and reaches Cd in M. As a result, from M starts the 
backward shock Cm given by 

Cm- x{t) = a{TZBit,x{t)) ,p{pi)) , x{tM) = xm- 

If tF is the solution of 

rCuit) rO 

F: I w{x) TZB{t,x) dx + p{pi) / w{x) dx = ^2, 

■J —iw JC]^i{t) 

t > tM , CM{t) > —iw , 

with Ciitp) < — itoi then the fall in the efficiency of the exit to p2 affects the flow 
and from F starts a shock Cf with constant speed a {p{pi) -, p{p2)) < that reaches 
Cm in N . From A'^ then starts the backward shock Cn given by 

Cjv: x{t) = a(TZB{t,x{t)) ,p{p2)) , x{tN)=XN- 

We assume that Cn and Cl meet in O with xq < —iw- Then from O starts a 
forward shock Co- Observe that ^(t) = p{p2) for any time t between tp and the 
time at which Cq crosses x = — i^, see Fig. [8] (b). After that, the map [t 1-^ C{t)] 
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starts to decrease and, consequently, the efficiency of the exit increases at time to 
when ^(tc) = ^2 and then again at time tfj when ^{tn) = Ci- Observe that the 
shock Co moves faster after its interaction with the two rarefactions started from 
G and H and that it finally reaches x = at time tj, that corresponds to the 
evacuation time. 




0,2 0.4 0.6 0.8 1.0 20 40 60 80 " PL 

(a) The functions [p >-^ f{p)] (b) The function [t h->- 5(t)]. (c) In gray the region C and in 
and [£. 1^ p{0]- white the region A/". 



Figure 8. The above figures refer to Sec. O 

Fig. [7] corresponds to a linear weight function w{x) = 2i~^(i^ + x), a normalized 
fiux f{p) = p{l — p) (namely the maximal velocity and the maximal density are 
assumed to be equal to one) and to the values 

Po = 0.21, pi = 0.168, P2 = 0.021, ^1 ~ 0.566, 6 ~ 0.731, 

XA = —5.75, xb = -2, iw = 1, xo = -2, xm ~ -0.4002, 

tc = 2, tD = 5, tE ~ 9.651, to ~ 85.045, tj ~ 87.498. (32) 

In Fig. [8] (b), is represented the corresponding map [t i-> S,{t)]. In Fig. [8] (c), 
we represent the region C, in gray, and the region Af, in white, as introduced in 
Definition [2j Notice that for this choice of / and p the region C4 happens to be 
empty. 

6. Technical section 

Lemma 3. Consider the family of scalar conservation laws 

dtu + d^f''{u) = (t,x) G [0,r] x]-oo,0[ , (33) 

and assume that /" converges uniformly on compacts to f as n goes to infinity. Let 
p^ he a sequence of Kruzkov entropy weak solutions to ([33]) . If p^ converges a.e. to 
pm[0,T]xR, i/ien/"(p"(-,0-)) converges weakly to f {p{-,0-)) in'L^{[0,T];R). 

Proof. First of all notice that by straightforward passage to the limit p is a Kruzkov 
entropy weak solution to 

dtu + d:rfiu) = (t,x) e [0,r] x]-oo,0[ . (34) 

Further, />" and p admit strong left traces on [0,T] x {0}, see Refs. [26], [33] . 
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Now let 6e be as in Q. Choose any 9,{p £ C^ 
0, ip{0) = 1 and observe that 



^1-^+ 



) such that 6(0) = e{T) 



4>it,x) 



6e{z) dz 



t-T+e 



x+l/e 
x+e 



5e{z) dz 



e{t) ^(X) 



is a C^-map with support in [0,T] x [— 1/e, 0] and, as e goes to zero 

(/)(0, x) = = (t){T, x) , dt4>{t, x) ^X[Q^ 2^] X ]_oo^ q] *^*' ^^ ^'^^^ "^^^^ ' 

(/)(t, 0) = , aj.(/>(t, x) ^X[Q^ 2^] X ]-oo, 0] *^*' ^) ^*^*) '^'^^^ 

By hypothesis /j" is a weak solution of ()33p . therefore 

[p"ai0 + r(p")a,,^] dxdt = o, 





and letting e go to zero we have 

JR- '- 



dx dt 



r(p"(t,o-))0(t)dt. 



By definition p is the strong limit in Lj-^^ of the sequence p^ as n goes to infinity 
and f^ ^ f uniformly on compacts. As n goes to infinity, the left hand side of 
the above equation converges to 



pipe + f{p) 9 ^ 



Arguing as above, we also find that 

p^9 + f{p)e^ 







dx dt 



dx dt . 



fipit,o-)) eit) dt , 



because p is a weak solution to ([M]) . Therefore Jq f"'{p"'{t,0—)) 9{t) dt converges 

to /g f{p{t,0—)) 6{t) dt as n goes to infinity, and the weak limit of /"'(/^"'(t, 0— )) 
equals /(p(t,0-)). D 

6.1. Proof of Proposition [31 For any £ G N and t G ]0, At/,], we know by Ref. [8] 
that the function /0^1]^(t) is piecewise constant with jumps along a finite number 
of polygonal lines. Therefore, for any t in a sufficiently small left neighborhood of 
i, we can write, by (|2ip 



p"''^(t + Mt,,x)=/>^;'\(t,x)= Yl P-U 



i&jl 



n,h 



n,h f 



n,h/,\ n,h , ( n,h n,h \ /. 7\ 

S^;Xt) = x_v + a [p^^i^pj^i^^j [t-t] , 



r x 



(35) 
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where J7"' C 



X 



h p, I n,h n,h\ r\ ^^h p, n.h /,\ , 'n,h/,\ 

= 0, fT(^p_'_o,p_'ij = 0, sJ_o = 0, s_;._i(t) < sj.(t), 

i)- 



/5- ■ G -M^ />_',■ / P_;-+i for any . / and f{p''_:l) = fip'^^) < p^{ 
Introduce the notation 



/-•• = r 



P-4 I ' 



* 



n,h 



^{P 



n,h 



P+ =P[P+ ,, 



~n,h - I n,h^ 

P+ =P[P-C ;, 



n,h 
p^ =p 



h f'^n,h 



m) 



ATy'^' (t) 



T^''^ (t+) 






it- 



Observe that by definition: 

if t < At/j, then 



n,h 



n,h 

p+ ; 



n,h 



n,h 



rn,h _ rn,h ^ n,h_ 
j-,0 — J~,l - P- ' 



n,h n,h 



rn,h 



-5-2-^ fip); 



n,h 



n,h 



(36) 
(37) 
(38) 



if pj 1 < p_'o, then /_'q = /_'i 

We have to distinguish the following two main cases: 

(U) either i= Ath, and in this case T'^{i+) - T'^it- 
(I) or t 7^ At/i, and in this case r^{i+) = T^(t—) and p"'™ = p"'"" by 

Let us consider more in detail the case t = At/j. If the efficiency of the exit 



n,h 



n.h 



does not change, i.e. p_' = p_^! , then by hypothesis H2, p"'' (t) is still given 
by ()35p for t lying in a sufficiently small right neighborhood of (l + l)Ath, and 
therefore AT^' (t) = —5 • 2~^ f{p)- If the efficiency of the exit changes, then 
2~^f{p) by Lemma [1] and we have to distinguish the following 



n.h 



P- 



n.h 
P+ - 

cases: 
(Ul) If the efficiency of the exit grows, then p"' = p^ + 2~^ f{p) and /"'q 
/"'i < P+ by ([37|) . There are two possibilities: 



n,h 



n,h 



(Ula) If /7^'q < p"^ I, then the solution does not change its expression, see Fig. |9l 
left and center, and 



^n,h 



AT"''^(t)=4 f{p)-p 



n,h 



/(p-)-P-'' 



5-2-^f{p) = -9-2-^f{p). 



n,h 



f^\ = p"' by ([38]) . For t > t sufficiently small. 



(uib) if/,:i';;</,:i';^,then/„o 

the solution contains a rarefaction between p"' q and p^ on the left of the 
constraint, a nonclassical shock between p^ and p^ at the constraint and a 
rarefaction between p^ and p"' -|^ on the right of the constraint, see Fig. [9l 
right. Therefore 



n,h. 



AT"'"(t) 



n,h 



n,h 
PJ 



f{p)-P 



+ 2 

n,h 



m-Pl' 



+ 



n,h n,h 

P+ -pJ 



5.2-'^/(p) = -5.2-V(p). 
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•n,.H 

P-,0 



77 ,h 

P-,0 



n .h 



iz ,h 





Figure 9. Interactions of the type (Ul). 



n,h n,h 



(U2) If the efficiency of the exit decreases, then p_,_' = pj —2 f{p). Differently 



n,h 



from the case (Ul), we do not know a priori whether p^ is less than 
/_'q = f_\ or not. Therefore there are three possibilities: 



71 ,h,^ 




n.k 
>'-,0 




Figure 10. Interactions of the type (U2a). 



(U2a) If /9^'q < p"^ ^ and /"'q = f^\ < p"' , then the solution does not change 
its expression, see Fig. [TUl and 



^n,h/ 



n,h 



AT^''^(t) = 4 /(p)-pr -4 /(P)-P-'" -5-2-^ f{p) = -2-'fip). 
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''-.0 




■n ,h 
''-.0 




P-,0 




Figure 11. Interactions of the type (U2b), (U2c). 



n,h 



ri,h 



n,h 



rn,h 



(U2b) If /9"'q < p"'^ and p^' < /"'q = /"'^, then after time t the solut 



ion 



n,h 



performs a shock between /9_'q and p^ on the left of the constraint, a non- 
classical shock between p^ and p^ at the constraint and a shock between 



n,h 



p^ and p^ I on the right of the constraint, see Fig. \TT\ left and center. 
Therefore, if /7^'q / /o"'ii then 



AT^''^(t) =2 


2/(p) - r;t - 


Pt 


+ 2 


/(p)-pf 


-2 


/(p) - /-; 


1 


-4 


'f{p)-p^\ 


-5-2-'^ /(p) = -2-V(/>) , 


while, if /9"'o = /9'!;^i, then 


AT^''^(t) = 


'2fip)-f-i-pt 


+ 2 


/(p)-pf 


+ 








— I 


i 


'f{-p)-p''- 


h" 





5-2-' 


^ /(p) = -2 


-V( 


O). 





(U2c) If p"'-^ < pi'o, then p"' < p"' = /"'q = /"'-^ by ([55|) . and after time 
t the solution performs a shock between p"'q and /5'|i on the left of the 
constraint, a nonclassical shock between /5'j: and p^ at the constraint and 



n,/i 



a shock between p]^ and pj^ on the right of the constraint, see Fig. [TTl 



CROWD DYNAMICS AND CONSERVATION LAWS WITH NON-LOCAL CONSTRAINTS 31 



ight. Therefore 












AT'^^ii) = 


'pn,h_^n,h 


+ 2 


/(P)-Pf 


+ 




-2 


/(/>) - 


n,h 


-5-2-^ f{p) = 



n,h n,h 

P- -P+ 



-2-V(p) 



In conclusion, for the case (U) we proved that AT^' ({£ + l)Ath) < —2 ^f{p) < 
-2-"/(^) for any £ G N. 

Now, assume that we have t G ]0, Atfi[. If at time i no interaction occurs, then 
AT^' (t) = 0. The remaining part of the proof consists in a detailed study of all 
possible interactions. We start with the most classical case when the interaction 
occurs away from x = 0. 

(10) If two waves, respectively between p^^^ and p"' and between p^' and p^^i, 
interact away from the constraint, then at least one of the two waves has 
to be a shock and, in any case, the resulting wave is a shock between p^^^ 



and p^^i- Therefore 



n,h/ 



AT^'"(t) 



^ 



n,h 



^ 



n,h 



■,i+l 



* 



n,h 



A-1 



* 



n,/i 



* 



n,h 



* 



n,h 



,i+l 



<0 



and the number of waves after the interaction diminishes. 
We now study the case in which a rarefaction reaches x = 0. 
(II) If a rarefaction reaches x = from the left, then /9^'q < p^ _i < p and 

/"'q+2~" = /"'_^. In particular, the solution cannot perform a nonclassical 

shock at X = before the interaction. There are therefore three possibilities: 

(Ila) If p"'q = p"'j^ and /"'_^ < p± , then the rarefaction crosses x = and 

(lib) it yO_ Q = p_ ]^ and p^ < J_ _i, then p^ = j_ q = j_^^ because n > n 
and for ()37p . After the interaction, the solution performs a shock between 
p"' _-^ and p^ on the left of the constraint and a nonclassical shock between 



-h 



n,h 



p^ and pj I at the constraint. Therefore 



-n,h 



AT^''^(t)= 2/(p)-/. 



n,/i 



n,h 
■P± 



+ 2 



n,h 



/(P)-P±' 



n,h 



f{p) -P± 



n,h 
P± 



f- 



n,h 



r'n,h 

'-Vip) ■ 



n,h 
■P± 



(lie) If p"'q 7^ p"'n then after the interaction the solution performs a shock 
between p"' __-|^ and p"' ^^ on the left of the constraint. Therefore 



n,h 



cn,h 



AT:^\t)= 2f{-p)-r-\-r-i 



-2 



f{-p) - /";i 



rn,h pn,h 



J-,1 ~ /-,-! 



)1— n 



/(Z')- 



(12) If a rarefaction reaches x = from the right, then the analysis turns out 
to be completely analogous to that for the case (II), and therefore it is 
omitted. 
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We now study the cases when a shock reaches x = 0. 
(13) If a shock reaches x = from the left, then p"'_x < mm{p^ q, p}, f"^'_i < 
/"'q = f^\ < p^ by ([37|) . There are three possibiHties: 
(I3a) If p"' Q = p"2 1, then the shock crosses x = and we have AT J' (t) = 0. 
(I3b) If /9"'q < p"j I, then after the interaction the solution performs a shock 



between /o"'_x and p^ ^ on the right of the constraint. Therefore 



^n,h 



n,h 



cn,h 



AT^'''(t) = 2f{p) - rj;[ - rj^ 



rn,h rn,h 



fip) - f-i 



. 



/"TO \ Tr ^i^ ^ ^M j-i ^M £n,h £n,h i 

(I3c) If pj 1 < p'q, then p±' = f\ = f_\ by 



, and after the interaction 



the solution performs a shock between p"' _^ and p"' ^^ on the right of the 
constraint. Therefore 



n,h/ 



AT"'''(t) 



n,h 
P± 



rn,h 
J- -I 



+ 4 



m-p^±' 



r, £/ -\ n,h £11, h 



/(P)-P±'' 



. 



(14) If a shock reaches x = from the right, then the analysis turns out to be 
completely analogous to that for the case (13), and therefore it is omitted. 
This concludes the proof of Proposition [3l 



6.2. Proof Proposition [4l We have to prove that the conditions given in Def- 
inition [1] . 
Consider 



inition [T] are satisfied. Fix cj) G C^(]R^;M+) and k € [0,R]. Let 5^ be as in 



t-iAt 



h 



Mt, x) = y] 4it - ^^th, x) = y] (t>{t, x) I 6e{y) dy 

^gN Jt-{e+l)Ath+e 

,), (/)f(-,x) has support in [0, At/^] for any 
X G M, and as e goes to zero 

<Pe{t,x) -^ <p{t,x) , 



and observe that 0£,</'e G C^(M2 



dt(t)e{t, x) -^ dt4>{t, x) + ^ (pit, x) 
dx(j)£{t,x) -^ dx<i){t,x) 



D 



^lAth (*) ~ \i+l)Ath (*) 



By construction, since each [t >-^ Pgliit)] is an entropy weak solution of (|19p in 
the sense of Definition [TJ we have 



n,h 



p^'^{t,x) — k dt(j)^{t,x) dx di 



At;,, 



+ 



+ 2 



n,h 



n,h 



sign(p^;'\ (t, x) -k) r[ p^;'\ (t, x) - r{k) d^cPi (t, x) dx dt 



Ath 



^n,h 



it) 



/"(p) 



rik)4{t,o)dt>o 
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By sumiiiing over ^ and letting e go to zero, by ()2ip we obtain that 

'0 

+ /" |/'m(0'^)-^ c^{i/^th,x)dx- I p";'\(At;„x)-A;,/)((^ + l)Atft,x)dx 

+ [ '[ sign(p^;'\(t, x) - k) \r (p"i(i, x)) - rik)] d,<P it + iAth, x) dx dt 
Jo Jr l \ / j 



/ / /9"l^i(t, x) - A; at</>(t + lAth, x) dx dt 
Jo Jm. 



+ 2 



At,, 






r{k)<p{t+iAth,o)dt 



p"^''{t,x)-k dt(t>{t,x) dxdt 



+ 



+ 2 



sign{p'''\t,x) - k) [r {p'''\t,x)') - rik)] d,cP{t,x) dx dt 



rip) 



r{k) (^(t,0) dt+ \p}^{x) - k\ 0(0, x) dx . 



Finally, by construction, /" (/)"•''(*, 0±)) < p^ {E''^^{t)) for a.e. t G M+, and this 
ends the proof of Proposition [H 

6.3. Proof of Proposition [6l We list here two basic properties which will be of 
great help in the following case by case analysis. 

First, by definition ([2]), ^(0) = pi and the map [t i— )• ^(t)] is continuous. Thus, 
by assumption (P2) we have that for any t > sufficiently small 

bpl: if e(t) < PL, thenj>(e(t)) =p(pl-); 
bp2: if PL < m, then p(e(t)) = p{pL+). 

The case ^(t) = pL is somehow special and has to be studied separately for each 
specific case. 

Second, when the solution is nonclassical, due to the finite speed of propagation 
of the waves, the assumption (P2) and properties bpl and bp2, we have 

npl: if n [pL,p{p)] (x) = PL for x < 0, then p = /(pl) G [p{pl+),p{pl-)]; 

np2: if p / /{pl) and pL < p{p)-, then p = p{pL+)] 

np3: if ^ 7^ /(pl) and p{p) < pL, then p = p{pL-); 

np4: if p is continuous in pL, namely p{pL—) = piPL+), then p = p{pl)- 

Now we start the description of the possible cases and we proceed as follow- 
ing. First, we show that for any initial datum satisfying (Ci), i = 1, ... ,5, the 
problem actually has a unique solution and that the solution is classical. Second, 
we take into consideration the corresponding case (Ni), for which we prove that 
the classical solution is not suitable and that there exists a unique nonclassical 
solution. 
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It is important to stress that in general the solutions to the constrained Riemann 
problem ()30p are not self-similar, see Example[TJ All the cases listed below describe 
self-similar solutions because we let the solutions to evolve only on a small interval 
of time. 

(CI) In this case [{t, x) i-^ T^Ipl, PR]{x/t)] performs a shock with negative speed 
'^{pl-.Pr) and satisfies ()30bp because f{pR) < p{pL+) and p{i{t)) = 
p{pl+) by bp2. Assume that there exists a nonclassical solution of the 
form ([3T]l . Observe that the assumptions pi < Pr and /(pr) < Kpl) 
together imply that p < pR. Then p{p) < p < Pr and TZ[p{p), pr] 
is given by a shock with non negative speed if and only if p < /(pr), 
or equivalently, pR < p{p). As a consequence, p < /(pr) < Kpl), 
PL < PR ^ p{p) and by np2 p coincides with p{pL+)- In conclusion 
we have p < /(pr) < p{pl+) = P, namely /(pr) = p and the nonclassical 
solution coincides with the classical one. 

(Nl) In this case [{t,x) t-^ TZ[pL, PR\{x/t)] does not satisfy ()30bp because 
Kpr) > PiPL+), see case (CI). Therefore, there does not exist any 
classical solution and we can consider only nonclassical solutions of the 
form (j3ip . If p is continuous in p^, then by np4 we have that p = p{pl)- 
If p experiences a jump at pi then, one may wonder which value in 
[p(pl+)iP(pl~)] has to be chosen as p. As in the case (CI), the as- 
sumptions imply that p < pR and then that p{p) < pR and p < /{pr)- 
Then p is strictly smaller than /(pl) and p{p) > pi. As a consequence, 
property np2 forces us to choose the unique possible value of p, which is 

P{PL + )- 

(C2) In this case [{t,x) i— )• TZ[pL, PR\{x/t)] performs a shock with non negative 
speed cr{pL,pR) and it satisfies (|30bp because /{pl) < p{pl+)- Assume 
that there exists a nonclassical solution of the form ()3ip . Observe that 
the assumptions pi < pR and /(pr) > Upl) together imply that p > pL- 
Then p{p) > p > pi and TZ[pL, p{p)] is given by a shock with non positive 
speed if and only if p < /(pl)- Thus p < /(pl) < p[pl+) and this implies 
by (|31cp that p = /(pl) = p{pl+) and that the nonclassical solution 
coincides with the classical one. 

(N2) In this case [{t^x) ^^ lZ[pL, PR\{x/t)\ does not satisfy (|30bp because 
Upl) > p{pl—), see case (C2). Therefore, there does not exist any 
classical solution and we can consider only nonclassical solutions of the 
form (f3T|) . As in the case (C2), the assumptions imply p{p) > p > p^. 
Furthermore, by (|31cp we have p < p{pL — ) < Upl), and as a consequence, 
property np2 forces us to choose p = p{pL+)- 

(C3) In this case [(f, x) i-> TZ[pL-,PR\{x/t)\ performs a possible null rarefaction on 
the right of the constraint and it satisfies ()30bp because /{pi) < PiPL+)- 
Assume that there exists a nonclassical solution of the form (|3ip . Since 
PL < P ^ p{p)i T^^PLi p{p)\ is given by a shock that has non positive 
speed if and only \i p < /(pl)- Therefore p < /(pl) < PiPL+) and 
this by (|31cp implies that p = /(pl) = p(pl+) and that the nonclassical 
solution coincides with the classical one. 
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(N3) In this case [{t,x) ^ 7l[pL,pB\{x/t)] does not satisfy (|30bp because 
Upl) > p{pl—), see case (C3). Therefore, there does not exist any 
classical solution and we can consider only nonclassical solutions of the 
form ([HT|) . By hypothesis and (|31c|) we have /(pl) > p{PL—) > P- There- 
fore PL ^ P < p{p) and by np2 we have p = p{pL+)- 

(C4) In this case [(i,x) H- TZ[pL-,PR\{x/t)\ performs a rarefaction with speeds 
between \{pl) < and X{pr) > and it satisfies (|30bp because f{p) = 
p{pL+) implies that p{p) = f{p) for all p < pi. Moreover, it implies 
also that p is continuous in pi and therefore, by np4, any nonclassical 
solution of the form ([3T]) must have p = p{pl) = /(p); but in this case the 
nonclassical solution coincides with the classical one. 

(N4) In this case [{t,x) H- TZ[pL, PR]{x/t)] does not satisfy (jSObp because f(p) > 
PipL—), see case (C4). Therefore, there does not exist any classical solu- 
tion and we can consider only nonclassical solutions of the form ()3ip . 
(N4a) By assumption and (jSlcp /(pl) < p{pl+) < P and therefore p{p) < 

PL and by np3 we have p = p{pL—)- 
(N4b) By assumption and (jSlcp f{pL) > p{PL—) > P and therefore p{p) > 
PL and by np2 we have p = p{pL+)- 

(C5) In this case [(t, x) i-)- TZ[pL, Pr\ {x/t)] performs a possible null rarefaction on 
the left of the constraint and it satisfies ()30bp because fipn) < p{pl—) and 
p{i{t)) = p{pl—) by bpl. Assume that there exists a nonclassical solution 
of the form ([3T]) . Since by assumption and ()31cp p > p{pL+) > fipi), we 
have p{p) < PL and by np3 p = p{pL—), but in this case the nonclassical 
solution coincides with the classical one. 
(N5a) In this case [{t,x) i-^ TZ[pL, PR\{x/t)] does not satisfy ()30bp because 
fiPR) > p{pl—), see case (C5). Therefore, there does not exist any 
classical solution and we can consider only nonclassical solutions of the 
form ([31]) . 
(N5b) By assumption and ()31cp . /(pl) < p{pL+) < P and therefore 

p{P) < PL and by np3 we have p = p{pL—)- 
(N5b) By assumption and (|31cp . f{pL) > p{pl—) > P and therefore 
p{P) < PL and by np2 we have p = p{pL+)- 



6.4. Proof of Proposition O 

(Rl) Any solution given by IZ* coincides on each side of the constraint with a 
solution given by the classical Riemann solver TZ. Therefore it satisfies the 
Rankine-Hugoniot jump condition along any of its discontinuities away from 
the constraint. Finally, by definition of p and p, it satisfies the Rankine- 
Hugoniot jump condition also along the constraint. 

(R2) It is clear by the proof of Proposition 

(R3) It is proved as in (Rl) since any classical solution is in BV. 

(R4) As was proved in Ref. [8], TZ* is continuous on C U A/". If {pl,Pr) is not in 
C L) J\f then p experiences a jump at S, = Pl- Therefore, the local in time 
solutions of the Riemann problem for the initial conditions (p/, + e, p/j) and 
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{pL — e, Pr) are different and only one of the two converges to TZ*[pliPr\ as 
e > goes to zero. 

(R5) We first stress once again that we can discuss the consistency property of 
our Riemann solvers only locally in time because, in general, the solutions 
may be not even self-similar globally in time. However, locally in time, the 
efficiency of the exit can be assumed to be constant and it is thus sufficient 
to proceed as in Ref. [H] . 

(R6) It is clear by the proof of Proposition O 

7. Further discussion on the model, conclusions and perspectives 

The present model does not take into account extremal cases. For instance, 
whenever a high density is approaching the exit the efficiency of the exit can 
become very small. As a consequence, even a small density of pedestrians may 
form a queue provided a sufficiently high density is approaching from behind. 
However, at least in this case, the low efficiency of the exit has not a "big" effect 
on the flow at the exit which is in fact "small" . Further investigation and modeling 
may be needed in order to deal with such singular effects. 

The planned forthcoming papers of the authors aim to generalize the present 
model to the initial-boundary value problem with non-local constraint, to code 
the resulting model, to simulate realistic evacuations, to state and solve optimal 
management problems, to reproduce the so-called Braess' paradox for pedestrian 
flows (e.g. the situation where a wider door placed before the exit door makes the 
evacuation faster), see Refs. [9], [TO], and to introduce further features of "panic" 
behavior, see Refs. [Ij, [13]. 
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